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ABSTRACT 

Context. The dust is an important tracer of dense interstellar clouds but its properties are expected to undergo changes affecting 
the scattering and emitting properties of the grains. With recent Herschel observations, the northern filament of the Corona Australis 
cloud has now been mapped in a number of bands from 1.2 iim to 870 iim. The data set provides a good starting point for the study of 
the cloud over several orders of magnitude in density. 

Aims. We wish to examine the differences of the column density distributions derived from dust extinction, scattering, and emission, 
and to determine to what extent the observations are consistent with the standard dust models. 

Methods. From Herschel data, we calculate the column density distribution that is compared to the corresponding data derived in 
the near-infrared regime from the reddening of the background stars, and from the surface brightness attributed to light scattering. We 
construct three-dimensional radiative transfer models to describe the emission and the scattering. 

Results. The scattered light traces low column densities of A v ~l m better than the dust emission, remaining useful to A v ~ 10 - 15 m . 
Based on the models, the extinction and the level of dust emission are surprisingly consistent with a sub-millimetre dust emissivity 
typical of diffuse medium. However, the intensity of the scattered light is very low at the centre of the densest clump and this cannot 
be explained without a very low grain albedo. Both the scattered light and dust emission indicate an anisotropic radiation field. The 
modelling of the dust emission suggests that the radiation field intensity is at least three times the value of the normal interstellar 
radiation field. 

Conclusions. The inter-comparison between the extinction, light scattering, and dust emission provides very stringent constraints on 
the cloud structure, the illuminating radiation field, and the grain properties. 

Key words. ISM: clouds - Infrared: ISM - Radiative transfer - Submillimeter: ISM 
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O ; The interstellar clouds are hierarchical structures where gray- ^^Ttetio^d S^ha^b^nSh. numer 

^ . nationally bound prestellar cores are found on the smallest Qus cos mological simulation s, and in other theoretical stud- 

. . scales. The study of the dense clouds is largely motivated by igs ( Burkert & Hartmann 2004) The filaments should f 

> this connection to star formation The filamentary nature of ment as dictated b me ^ Je an S condition, a proce ss ad- 

^ ■ interstellar clouds and the possible connection between fila- dressed by_many the0 retical studies (Inutsuka & Mivamall 1997b 

°»ent s and If^fg 10 " , Jb n^ L 1 ^° f 2?" g ^ 1 Evers 2009). In addition to the pure compressive instability, 

d i^ B c arnard llgUpssenkov 1952; Elme^^& Elmegreen| also the so _ called s ausage instability may sometimes play an 

[rMlSchneider&Ilmegreen 1979; Bally et al.||1987|). The re- aHHirinna1 ^ Fk L ra HZin l bm 9 l I m.t .Ln ,Ll I 



1, Introduction plosions and the radiation and stellar winds from massive 

stars. However, some filaments may also be formed directly 



cent ground-based studies and Hersc hel surveys have shown 



cloud filaments in exqu i site detail (Miville-Deschen es et al 
I201QI: I Andre et all l201Qt lArzoumanian et al. 201 1] iHill et al 



addi tional role (see iFischera & Martini I2012I iMcLeman et al.l 
I2012I) . Although the morphology of all the star-form ing clouds 
is not predominantly filamentary (lJuvela et al.ll2012l) . a global 
picture of star formation is emerging, where the turbulence 



I2011J Uuvelaet aL l |2010) and have demonstrated that pre- creates filamentS; the filaments become gravitationa i ly unsta . 

stellar cores and protostars are p r eferentially located along w and subsequentl fragment forming the cores that ma 

these sJ^u^J^M Emmm EE 'shchikovetalj ioIg ^ fed ^ ^ ^ ^ fiiaments ^ 

| JuvelaetaL | 2012|) Fila ments are a natural otone of inter- mass ^ can bg found ^ individual filaments while hi h 

stellar turbulence (e.g. | Padoan & Nordlun J M), with fur- mass §tars ^ ferentiall born at the intersections of sev- 

ther contnbutions from immediate triggering by supernova ex- eral filaments This scenario is supported by the observations 



Herschel is an ESA space observatoty with science instruments "22^ within the Herschel Gould Belt survey (e.g. |Andre etal] 



nerscnei is an m/\ space ooservaiory witn science instruments . "\ , , -. n . ■ 

provided by European-led Principal Investigator consortia and with im- I2Q10|; IKonyves et ai pOlfllMen'shchikov et al. 2010) and other 
portant participation from NASA. Herschel programs (IMolinari et al.ll2010l ISchneider et al.ll2010l : 
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Ngu venLuong et al.1 120111: iHill et al.1 120111: ISchneider et all 
2012; Juvela et al.l l2012h as well as by numerical sim- 



ulatio ns (Padoan & Nordlund 2011; Vazqu ezes emadeni et al.l 
1201 lblKlessenll201 ll i lBonnell et alJl201 lh ~ 

To determine the initial conditions for the star formation to 
occur, we need to measure the physical properties of the clouds, 
the filaments, and the cores. The mass distributions of the clouds 
can be measured with a number of methods. This is fortunate be- 
cause each method suffers from different sources of uncertainty. 

The dense clouds have been traditionally mapped in molec- 
ular lines. Line observations are invaluable because of the infor- 
mation they give o n the physical sta te, kinematics, and chem- 
istry of the clouds (iBergin & Tafallall2007l) . On the other hand, 
they do not always provide reliable or consistent estimates of the 
mass distribution. The line intensity depends on the local phys- 
ical conditions, mainly density and temperature, and the chem- 
ical abundances. If the lines are not optically thin, the radiative 
transfer effects introduce additional uncertainty. 

For the above reasons - and because of the development in 
detector technology and the appearance of new ground-based, 
balloon-borne, and space-borne facilities - the thermal dust 
emission has become important as a tracer of the densest clouds. 
Observations of the dust at submm wavelengths are of partic- 
ular interest because they are sensitive to the emission of the 
cold dust that may have been missed in earlier far-infrared stud- 
ies. With the knowledge of dust temperature and dust opac- 
ity, the column density can be estimated. However, the colour 
temperature obtained from observations is known to be a bi- 
ased e s timator of the average dust temperature dShetty et al.l 
I2009bt lJuvela & Ysardl l2012bh and additionally the value of 
the dust opacity is uncertain. Furthermore, there are clear in- 
dications that the dust opacity is not constant, the variations 
being a likely c onsequence of grain coagulation and aggre- 
gation processes (Cambresv et al. 2001; del Bur go et al.l 12003: 
iKramer et all 120031: iLehtinen et al.1 120071) . Similarly, the spec- 
tral index appears to vary from region to region and this in- 
creases the uncertainty of the colour temperature and col- 
umn density estimates. The observations point to a negative 
correlation between the colour temperature T& and the ob- 



served spectral index (IDupac et al. 2003; Desert et al. 2008; 
Anderso neTaflkOlOl: I Paradis et al. 2010; Venezia nietaDl2010 : 
Planck Colla boration et al.ll201 lbtlArab et al.ll2012h . Such varia- 
tions can affect the accuracy of the column density estimates de- 
rived from the dust emission. However, the evaluation of errors is 
complicated by the observa tional noise ( also through its effect on 
the T d - ytfobs relation, see IShettv etaD (l2009ah : ljuvela & Ysardl 



2012al) and the un known line-of- sight tem perature variations 
(Shett y etal]l2009bt I Juvela & Ysardl l2012al) . One of the main 
advantages of submm observations of dust emission is the large 
range of column d ensities probed. Wi th the Herschel satellite 
SPIRE instrument (Gri ffin etall 1201 Oh one can probe regions 
from Ay ~ l m to Ay ~ 100 m with spatial resolution of ~ 30" . 
However, the accuracy will depend on th e line-of- sight tempe r- 
ature variations and other opacity effects (Mali nen et al.ll201 H) . 

Complementary information can be obtained from near- 
infrared (NIR) observations, from reddening of the light of 
background stars, and from the measurement of the scattered 
light. The observed near-infrared colour excesses can be con- 
verted to column density or to provide estimates of the visual 
extinction Ay using tech niques such as the NICER algorithm 
(lLombardi & Alves 2001). The calculations are based on the as- 
sumption that the average intrinsic colour of the stars is constant 
and, therefore, reliable estimates are obtained only as an average 
over a number of stars. Therefore, the extinction must be esti- 



mated for larger pixels, as an average over -10 or more back- 
ground stars. This reduces the resolution of the extinction maps 
and can cause bias in the presence of strong Ay gradients. With 
dedicated observations it is possible to achieve a spatial reso- 
lution of some tens of arc seconds, which is comparable to the 
resolution of the submm surveys. The range of Ay values probed 
extends from Ay ~ l m to ~ 40 m , depending on the stellar density 
and the depth of the observations. 

If the NIR observations are made using an ON-OFF mode 
the diffuse signal is preserved and it is possible to measure the 
light scattered by the same dust particles that are responsible for 
the NIR extinction. If the cloud is optically thin at the observed 
wavelengths and if the cloud does not include local radiation 
sources, the signal will be proportional to the column density 
(lJuvela et al.1 12006). The main advantage of the surface bright- 
ness measurement is the high, potentially ~ 1" spatial resolu- 
tion. Even in normal NIR observations at mid-Galactic latitudes, 
the scattered light will provide better resolution than similar ob- 
servations of the background stars. The method remains useful 
in the range Ay = 1-1 5 m . At higher column densities the signal 
saturates and the surface brightness will depend on the structure 
of the cloud. With multiwavelength measurements (e.g., J, H, 
and K bands) it becomes possible to use the intensity ratios to 
measure and to partially correct for the saturation of the signal 
(lJuvela et al.ll20 06. 2008, hereafter Paper I and Paper II, respec- 
tively). The observed surface brightness depends on the grain 
properties and upon the intensity of the radiation field at the 
location of the cloud. Therefore, combining surface brightness 
data with other column density tracers one can also provide con- 
straints to these parameters. 

In this paper we will examine a filament in the northern 
end of the Corona Australis molecular cloud. The filament con- 
tains a dense clump with a central Ay previously estimated to 
be at least ~ 30 m . The source has been studied with the help 
of NI R observations us i ng b ackground stars and the scattered 
light (lJuvela et al.ll2006l 120081 120091) . The region has now been 
mapped with the Herschel satellite and these data provide the 
opportunity to measure the structure of the central part of the fil- 
ament where the column density was too high to enable reliable 
estimates either with the NICER or with the surface brightness 
(scattered light) methods. This provides an excellent opportunity 
to compare the performance of the three column density tracers 
and to use the comparison to look for indications for changes in 
the dust properties. 

The structure of the paper is the following: In Sect. [2] we 
describe the observations and the main properties of the data. 
In Sect. [3] we present the direct results derived from each data 
set and make the first comparison between the column density 
estimates. We continue by constructing radiative transfer models 
for the NIR data and the sub-millimetre data separately. These 
calculations are described and compared in Sect. [4] Our final 
conclusions are given in Sect. [5] 



2. The observations 

2.1. Near-infrared data 

The NIR observations were made with the SOFI instrument on 
the NTT telescope, Chile. All observations were carried out in 
ON-OFF mode. The northern part, an area of 4'x8' was observed 
in 2006. The integration times on the ON-fields were 1.0 hour in 
the J band, 1 .5 hours in the H band, and 4.5 hours in the Ks band. 
The map was extended in June 2007 with two 4' x 4' fields that 
cover the southern side of the filament. The integration times 
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Fig. 1. A cross section of the northern Corona Australis filament according to the different tracers. Maps are shown for the visual 
extinction derived from the reddening of the background stars (frame a, 40" resolution), the 250yum surface brightness (frame b, 
18" resolution), the column density derived from the Herschel data (frame c, 40" resolution), and the intensity of the scattered light 
in the K s band (frame d, 4" resolution). For the plot the data have been adjusted so that the median signal is zero within the reference 
area located in the northern end of the maps (see text). The dashed rectangle in frame d denotes the area that will be used to measure 
the filament profile. 



were 1.5 hours for the Ks band and 0.5 hours in the J and H 
bands. On the basis of the 2MASS data dSkrutskie et al.ll2006l) . 
the OFF field was estimated to have an extin ction below Ay = 
0.5 m . For the details of the observations, see I Juvela et al.l (2008) 
an d lJuvela et al.l(l2009h . 

The observations were calibrated using the 2MASS stellar 
photometry and that calibration was carried over to the surface 
brightness. For the final surface brightness map the southern ON 
fields were mosaiced together with the northern part, this requir- 
ing a small additive adjustment. The observations do not provide 
an absolute measure of the surface brightness. Therefore, we 
present the results relative to the signal within a 1 ' diameter aper- 
ture centred at coordinates Of 2 ooo=19 h m 43.5 s , £ 2 ooo=-36°49'36" 
(see Fig.QJl). Based on the 2MAS S data, the visual extinction in 
this refere nce area is ~ 1.8 m (seeUuvela et al] |2009). 

In lJuvela et al.1 (l2Q08h the NICER method 
Combardi & Alves] l2001h was used to convert the mea- 
sured stellar colour excesses to estimates of extinction. This 
extinction map is reproduced in Fig. [T^ with Ay values relative 
to the above described reference area. The frames b and c show 
for the same area the 250 jim surface brightness and the column 
density derived from the Herschel data. 

2.2. Herschel observations 

The Corona Australis cloud h as been mapped a s part of the 
Herschel Gould Belt Survey (lAndre et al.l boiOh . The entire 
cloud is covered by parallel mode maps that also cover the north- 
ern fila ment. The observatio n were made with the PACS instru- 
ment (Po glitsch et al.|[20Toh at wavelengths 70/mi and 160/mi 



1 PACS has been developed by a consortium of institutes led by 
MPE (Germany) and including UVIE (Austria); KU Leuven, CSL, 
IMEC (Belgium); CEA, LAM (France); MPIA (Germany); INAF- 
IFSI/OAA/OAP/OAT, LENS, SISSA (Italy); IAC (Spain). This devel- 



and with the SPIR^l instrument (iGriffin et al.ll201Qh at wave- 
lengths 250yum, 350yum, and 500 yum. The observation consist of 
two orthogonal scans with observation identifiers 1342206677 
and 1342206678. We use data provided by the Herschel Gould 
Belt Consortium. We use the SPIRE maps produ ced with the 
naive mapping routine in HIPE software ( Ott 2010). The 160 yum 
PACS map is made using Scanamorphos v. 16 (Roussel 2011, 
submitted), starting from data processed to level 1 in HIPE. The 
70//m data were not used, since at that wavelength the emis- 
sion is not produced by the large dust grains that are responsi- 
ble for the longer wavelength emission. For the Herschel data 
the zero point of the intensity scale was estimate d by compari- 
son w ith data from the Planck satellite survey (see lBernard et al.l 
l2010h . These zero points are used in the derivation of the col- 
umn density estimates. However, in the plots (e.g., FigQJ)) and 
in comparison with the other tracers, the zero level is set using 
the reference area described in Sect. 12.11 

2.3. Laboca 870/um data 

The Corona Australis filament was mapped in 2008 with 
the APEX/LAB OCA instrument. This instrument operates at 
a wavelength of 870 yum, and has the beam size ~ 19". 



opment has been supported by the funding agencies BMVIT (Austria), 
ESA-PRODEX (Belgium), CEA/CNES (France), DLR (Germany), 
ASI/INAF (Italy), and CICYT/MCYT (Spain). 

2 SPIRE has been developed by a consortium of institutes led by 
Cardiff University (UK) and including Univ. Lethbridge (Canada); 
NAOC (China); CEA, LAM (France); IFSI, Univ. Padua (Italy); IAC 
(Spain); Stockholm Observatory (Sweden); Imperial College London, 
RAL, UCL-MSSL, UKATC, Univ. Sussex (UK); and Caltech, JPL, 
NHSC, Univ. Colorado (USA). This development has been supported 
by national funding agencies: CSA (Canada); NAOC (China); CEA, 
CNES, CNRS (France); ASI (Italy); MCINN (Spain); SNSB (Sweden); 
STFC (UK); and NASA (USA) 
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Observations consist of ~ 30' long scans that were made per- 
pendicular to the filament, and at up to 30 degree angles relative 
to the normal of the filament. The observations were reduced us- 
ing the Boa program, version 1.11 and the calibration checked 
with planet observations and with the help of the source S CrA 
that is located immediate south-east of the area covered by the 
NIR data. The final noise level in the central part of the map is 
bet ter than 20 mJy/be am. For details of the Laboca observations, 
see I Juvela etaD (l2009h . 

2.4. Mid-infrared data 

The area was mapped in Spitzer guaranteed time programs with 
both the IRAC instrument (3.6, 4.5, 5.8, and 8.0/im) and the 
MIPS instrument (24, 70, and 160/im). The MIPS data were al- 
ready compared to the LABOCA and NIR data in lJuvela et al.l 
( 2009). The images of the four IRAC bands are shown in Figj2] 
The 5.8 and 8.0/mi surface brightness is affected by the nearby 
stars and does not trace the column density. The 3.6 and 4.5 jim 
data will be used in Sect. 13.31 check for signs of enhanced mid- 
infrared scattering. 

We will use also data from the WIS^l survey that have re- 
cently become public. The WISE data include s the wavelengths 
of 3.4,4.6, 12, and 22//m dWright et al. 11201 Oh . 

3. The results 

3.1. The column densities 

3.1 .1 . Column density derived from the Herschel data 

The Herschel measurements of submm dust emission were con- 
verted to column density using the usual procedure that assumes 
a constant dust temperature along the line of sight and a constant 
dust opacity everywhere within the mapped area. For each map 
pixel, the observed intensities I v at 160//m, 250/mi, 350 /zm, and 
500 yum were fitted with a modified black body curve B v (Tq)^ 
to determine the colour temperature Tq. With knowledge of the 
dust opacity k (and the spectral index J3) the fit determines the 
column density. The column density of the molecular hydrogen 
can be written formally as: 

AW = , (1) 

where T is the assumed temperature, here equal to the colour 
temperature 7c, \i is the average particle mass per hydrogen 
molecule, and the dust opacity k is given relative to the gas mass. 

In the analysis, the surface brightness maps are convolved 
to a common resolution of 40" resolution before the deter- 
mination of Tq and N(H.2). The spectra are modelled assum- 
ing a constant value of the spectral index, /3 = 2.0 and a 
dust opacity k =0. 1 cm 2 /g ( v/1000 GKzf that is applicable to 
high density environments (Hildebrand 19831 iBeckwith et al.l 
1990). The same opacity law has been adopted in several pa- 
pers on Herschel results ([An dre et al. 2010; Konvves et al. 2010; 
Arzo umania n et al.l 120 1 it iPlanck C ollabor ation et al.l 120 1 1 a ; 
Juv ela et al.l l2012). However, for any particular sources, the val- 
ues of these parameters are not known to high precision and both 
the J3 and the k values introduce an uncertainty in the column 
density that may be several tens of per cent. 

3 The Wide-field Infrared Survey Explorer is a joint project of 
the University of California, Los Angeles, and the Jet Propulsion 
Laboratory/California Institute of Technology, funded by the National 
Aeronautics and Space Administration. 
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Fig. 3. The map of the dust colour temperature Tq. 

The calculated column density map is shown in Fig.Q]: and 
the colour temperature map is shown separately in Fig. [3] The 
minimum temperature at the location of the densest clump is 
Tq ~ 12.5 K. This value is obtained using the total intensity (i.e., 
the original zero point of the Herschel maps) without subtracting 
the signal in the local reference area (see Sect. 12.21 1. 

3.1 .2. The A v maps from NIR data 

The extinction derived from the colour excesses of the 
background stars was c alculated using the NICER method 
(Lombard! & Alves 2001). The resulting maps of visual extinc- 
tion, Ay ICER have a resolution of 30". A si milar map wi t h a 20 " 
resolution has already been presented by lJuvela et all (120081) . 
Here we settle for a more reliable but a slightly lower initial res- 
olution of 30" because the data will later be compared to lower 
resolution column density maps derived from the Herschel data. 
The extinction map covers the whole area included in the NIR 
observations (see FigOJ). The maximum extinction values are 
~33 m but because no background stars are visible at the centre 
of the dense clumps within the filament, this is only a lower limit. 
Furthermore, because of the strong Ay gradients th e values near 
the fi lament centre may be biased downwards (see lJuvela et all 

120081) 

In lJuvela et al] (120081) and lJuvela et al.1 (l2009h the surface 
brightness, assumed to consist only of scattered light, was con- 
verted to column density assuming analytical formulas I v ~ 
a v (l - exp(-b v Ay)) for the relation between the surface bright- 
ness and the extinction. This resulted in an extinction map A y CA 
that was consistent with the Ay ICER values up to ~ 10 m but gave 
higher values in the range of 10-1 5 m . It was not clear which of 
the two extinction estimates was more robust although simula- 
tions indicated that most of the difference could be explained by 
bias in Ay ICER . When the extinction approaches 20 m the surface 
brightness is strongly saturated and the simple analysis method 
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Fig. 2. Spitzer IRAC maps at 3.6, 4.5, 5.8, and 8.0 jim. The dashed line shows the outline of the regions covered by NIR observations. 
The brightest star east of the area covered by NIR observations is S CrA. The other bright source in the southern end of the images 
(at 6 2 ooo = -37°3'42") corresponds to the blended image of the stars HD 176269 and HD 176270. 




Fig. 4. A comparison of the Ay estimates derived from the sub-millimetre emission (frame a, the column density converted to Ay 
with the assumptions listed in the text), the reddening of the background stars (frame b), and the dust scattering (frame c). All maps 
have been convolved to a common resolution of 40". The frame d shows the NICER Ay estimates (black points) and the estimates 
A y CA derived from the NIR surface brightness as the function of the values from sub-millimetre data. The blue and the red points 
are, respectively, the A y CA values north and south of declination -36°55 / . The solid lines are least squares fits to the respective data 
points. The dashed line indicates the one-to-one relation. 



adopted in I Juvela et al.l (120081 120091) leads to diverging values. 
The saturation of the surface brightness and the dip at the centre 
of the dense clump are evident even in the K s band (see FigQ]i) 
where, compared to the V band the optical depth is lower by 
close to a factor of ten. This alon e confirms t hat th e central Ay 
must be higher than 20 m . As in I Juvela et aD (l2008h . the spatial 
resolution of this initial A y CA map is 10". 



3.1 .3. The comparison of the column density maps 

In Fig. |H we compare the column density estimates from the 
submm emission, t he reddening of the backg round stars (the 



NICER method by Lombard i & Alves (2001)), and from the 
NIR surface brightness data (lJuvela et al.l l2009). 



For easier comparison, the Herschel estimates of A/^Fb) were 
also converted to units of visual extinction, A y b_mm . The conver- 
sion was done using the canonical relation N(¥l2)— 0.94 x 10 Ay 
(IBohlin et al . 1978). That relation was determined for more dif- 
fuse lines of sight and may not be representative of the present 
field, but does at least provide a convenient point of reference. 
The map of A y b_mm (a scaled version of Fig. [T]:) is shown in 
Fig. [4^. The map of A y ICER was convolved to the resolution 
of 40". We have masked the areas near the borders of the ob- 
served area where the result of the convolution is not well de- 
fined (where more than 10% of the convolving beam falls out- 
side observed region). The A y ICER map is shown in Fig. [4J3. In 
the case of the NIR surface brightness, A y CA values remain un- 
defined in the central part of the filament where the simple an- 
alytical solution diverges. In the area where the A y CA estimates 
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exceed 25 m , before making the convolution, the data were re- 
placed with values read from the NICER map. After the con- 
volution, those areas are masked from the subsequent analysis 
(Fig. Ht). In the plots Fig. |4^-c the map zero points were set 
using the common reference region (see 12. It . 

The three Ay maps are morphologically very similar at low 
and intermediate extinctions. In the lowest Ay regions the abso- 
lute value of the visual extinction is of order of 1-3 mag . At those 
diffuse lines-of- sight the correspondence is best between the two 
NIR-derived maps because one is close to the sensitivity of the 
sub-millimetre data (or because the sub-millimetre maps may be 
affected by minor baseline uncertainties). 

Figure |4ji confirms the good correspondence between the 
tracers. Below Ay b_mm = 10 m the correlation coefficients 
coefficient is 0.96 for the A^ ICER -A^ h ~ mm relation with a 
least squares line A* ICER =-0.1 l + 1.76xA^ b_mm . For the relation 
Ay CA -Ay b_mm the correlation coefficient is r =0.86 but the pres- 
ence of two different relations is evident. In the northern half of 
the map (blue points in Fig.[4jl) the least squares fit gives a rela- 
tion A v CA =0.15+1.05xAy b ~ mm with a correlation coefficient of 
r =0.95. The slope s are consistent with the results presented by 
lJuvela et al.l (120081) where the A s y CA was below the A* ICER values 
up to ~ 10 m while the situation was reversed at higher Ay. On 
the southern side the relation is A V CA =0.81 +2.03 xA^ ub - mm with 
an even higher correlation coefficient, r=0.98. The main differ- 
ence compared to the northern side is the much steeper slope of 
the relation. The northern and the southern fields were calibrated 
separately, but a factor of two difference in the calibration is very 
improbable. Furthermore, the surface brightness values decrease 
to zero both at the southern and northern edges of the maps and 
the surface brightness values were matched at the centre when 
the fields were mosaiced together. Therefore, there is no room 
for an over 60% multiplicative error in the relative calibration. 
The difference is an indication of a change in the local NIR ra- 
diation field. Assuming that the grain properties are the same, 
and ignoring the fact that the scattering is mostly in the forward 
direction, this would suggest that the radiation field intensity is 
at least 50% higher on the southern side of the fi lament. This 
is qu alitatively consistent with the conclusions of I Juvela et al J 
(2009) that were based on the morphology of the NIR and the 
Spitzer dust emission maps. 

The average ratios A S Y CA /A™ h ~ mm and A% ICER /A™ h - mm are 
~ 1.5 or higher. The interpretation of A V CA is affected by the 
uncertainty of the intensity of the interstellar radiation field 
(ISRF) but that does not influence Ay ICE R It would thus ap- 
pear that ISRF is not very different from iMathis et all (Il983l) 
values north of the filament but are clearly elevated on the 
southern side. The Ay b_mm values appear to be too low by al- 
most a factor of two. Part of the difference could be caused 
the assumed scaling between N(H2) and Ay, but it is clear 
that especially the value of the dust opacity k has a large un- 
certainty. The Ay b ~ mm values were derived with an opacity of 
k =0.1 cm 2 /g (v/l OOOGHzy 3 that is higher than the value found 
in diffuse regions. iBoulanger et al.l (1 19961) derived for high lati- 
tude clouds a value of /c(850 yum)=0.005 cm 2 g" 1 which is -40% 
of the value given by the previous formula. Thus, with the diffuse 
medium k the Ay b_mm values would rise above the Ay ICER val- 
ues. This suggests that the actual value of k should be somewhere 
in between. However, it is also known that the colour tempera- 
ture used in Eq.[T] overestimates the mass averaged dust temper- 
ature and may lead to the underestimation of the c olumn den- 
sities (see IShettv et alJl2009bt iMalinen et aDl201ll: lYsard et al.l 
120121: lJuvela & Ysa rd 2012t|. If this effect can be quantified, it 



becomes possible to increase the lower limit the dust opacity. 
This is possible only with modelling and we will return to this 
question in Sect. HI 

3.2. The NIR and sub-millimetre profiles of the filament 

We measure the profile of the filament using the three tracers 
available to us. The sub-millimetre emission is the only one that, 
at least in principle, is capable of probing the column density 
distribution across the whole filament. Good Ay ICER estimates 
are missing for a small central part where no background stars 
are visible (roughly the white area in Fig. [4^). The NIR surface 
brightness data are reliable only in regions of less than Ay ~ 
10 m . Because we want to compare all the tracers, we concentrate 
on the narrow region that is marked with dashed lines in Fig. [TJi. 
We calculate the one-dimensional profile along the main axis of 
that region, averaging the data in the perpendicular direction. 

Figure [5] shows the resulting column density profiles that, 
like in Fig. HI are converted to units of visual extinction. The 
first frames show the data at the original resolution (cf Fig. [4]). 
There is an offset between the maps. In the northern end the 
^sub-mm decreases below 0.5 m while the other maps are, based 
on the large scale 2MASS extinction maps, were set at ~ 2 m . 
This is not likely to be caused by an uncertainty in the zero point 
of the Herschel surface brightness data. The missing ~ 1.5 m 
corresponds to a 250 /mi surface brightness of 28MJy sr" 1 , this 
assuming a temperature of = 17 K and a dust opacity of 
k = 0.005 x W850yum)- 2 cm 2 g- 1 (see Sect. 15X51 . Such a 
surface brightness offset would be comparable to the minimum 
signal found in the northern part of examined area. Figure 
again shows the discrepancy in the A y CA estimates. On the north- 
ern side these follow closely the shape of the Ay ICER curve but 
they are significantly higher in the southern end. In the plot 
we have not included the A V CA values for the central filament, 
Ay CA > 10 m , because of their large uncertainty near the region 
where the surface brightness saturates. 

The frame b of Fig. [5] shows the same data convolved to a 
common resolution of 40". The value at the offset of 10' has 
been subtracted, and the Ay b_mm values have been scaled by a 
factor of 1.3. With this scaling, the column density estimates 
are very consistent on the northern slope. For Ay ICER and the 
scaled version of Ay b_mm the match is good up to Ay = 15 m . 
However, the peak of Ay ICER is almost one minute of arc south 
of the Ay b_mm maximum. The Ay b_mm profile is strongly skewed 
towards the north and rough agreement is found again in the 
south end of the map. If we trust the NICER map, both A V CA and 
Ay b_mm show deviations on the south side of the filament and 
qualitatively both features could be affected by the asymmetry of 
the radiation field. Figure [5^ shows as dotted lines the minimum 
and maximum Ay ICER along the stripe for which the other curves 
show the average value. These indicate that there is a strong Ay 
gradient perpendicular to the stripe and parallel to the filament, 
especially close to the high column density clump. T his could 
cause significant bias in Ay ICER (see lJuvela et al.ll2008l) but can- 
not explain all the differences between Ay ICER and Ay b_mm , nei- 
ther in the levels nor in the profiles. 

3.3. Mid-infrared profiles 

Although we will not attempt to model the mid-infrared data of 
the Corona Australis cloud (Sect. [4]) it is also useful to examine 
the Spitzer IRAC and the WISE satellite observations of the fil- 
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Fig. 5. The column density profiles for the region marked with 
a dashed rectangle in Fig. [TJi. The left hand frame shows the 

A sub-mm ( tMck black lme)? A NICER (mm blue ^ A SCA 

(dashed red line, excluding the central filament where the NIR 
emission becomes partly saturated) at the original resolution 
(40", 30", and 10", respectively) without scaling. The x-axis is 
the offset along the stripe, from south to north. The dotted lines 
are the profiles of the maximum and minimum Ay ICER along the 
stripe. The right hand frame shows the same data at a common 
40" resolution, with values set to zero at offset 10' (additive cor- 
rection). The Ay b ~ mm values have been scaled by a factor 1.3 to 
match the other curves on the northern side of the filament. 




Offset (arcmin) 

Fig. 6. The mid-infrared profiles of the Corona Australis fila- 
ment. The offsets are the same as in Fig. [5] In addition to the 
Spitzer 3.6/mi and 4.5 jim data (blue and lines, respectively, and 
the right hand y-axis) the figure includes the profiles for the scat- 
tered light in the H and Ks bands (cyan and magenta lines, left 
hand axis) and the Ay derived from the Herschel observations 
(black dashed line, left axis). 



ament. Figure [6] compares the 3.6yum and 4.5 yum Spitzer IRAC 
profiles with the scattered light and the extinction derived from 
the Herschel data. 

The mid-infrared profiles show some resemblance to the NIR 
bands but follow more closely the column density profile of the 
filament as derived from the sub-millimetre observations. In par- 
ticular, the 4.5 yum emission peak coincides with the location of 
the column density maximum derived from the Herschel data. 
Compared to the estimated column density, both the 3.6 yum and 
4.5 yum data show significantly higher intensity levels on the 
south side. We compared the Spitzer and WISE data at the wave- 
lengths of 3.4yum and 4.6 yum. The WISE measurements are in 
complete agreement with the Spitzer data at the corresponding 
wavelengths except for the 3.4 /mi intensity that is slightly lower 
in WISE at the location column density peak (~0.18MJy/sr vs. 
0.21 MJy/sr) so that the intensity remains almost flat between 
offsets 3.0 / -5.5 / . Further outside the filament, the Spitzer ratio 
3.6 yum/4.5 yum and the WISE ratio 3.4 yum/4.6 yum show some in- 
crease. This could be caused by the effect the optically thick fil- 
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ament has on the local radiation field but the quantification of 
these effects would require separate modelling. 

The mid-infrared data have revealed in many dense cloud 
cores the presence of the "coreshine" phenomenon where the 
emission at -3. 5 um becomes brighter relative to the emis- 
sion at -4.5 /mi dSteinacker et al.ll2010l:IPagani et al.ll2010h . The 
changes are interpreted as a sign of an increase of dust grain sizes 
that leads to enhanced light scattering at wavelengths beyond 
3yum. The first detection of coreshine was made in the cloud 
LDN 183 using Spitzer data dSteinacker et al.ll2010b . Further de- 
tections have subsequently been made with WISE observations 
(I Juvela etal.1 12012b but LDN 183 remains the best example of 
the phenomenon. Figure [6] shows significant increase of the 
3.6yum/4.5yum ratio but only south of the column density peak. 
This is thus probably caused by a change of the radiation field 
rather than by grain growth. Regarding the overall asymmetry 
with respect to the filament centre, the 3.5 yum profile is rather 
similar to that of the K band. 



4. Radiative transfer models 

4.1. Model of the sub-millimetre emission 

We construct a three-dimensional model to explain the Herschel 
observations of 160 yum, 250yum, 350yum, 500 yum. In the present 
paper we examine only the main effects, including the pos- 
sibly anisotropic radiation field, using two dust mode ls. The 
first model corresponds to Ry=5.5 and is described in iDrainel 
d2003lfl . In the following, this dust m odel is called MWD 
As a p oint of comparison, we use the lOssenkopf & Henning 
(1994) dust model (in the following, the model OH) for co- 
agulated grains with thin ice mantles that have accreted in 
10 5 years at a density of 10 5 cm" 3 . At wavelengths below 
1/im we adopt the short wavelength extension discussed in 
IStamatellos & Whitworthl d2003h . The dust models differ with 
respect to the sub-millimetre spectral index f3 that, measured be- 
tween 250 yum and 500 yum is ~ 2. 1 for the MWD dust model and 
1.76 for the OH model. The optical depth ratios for extinction 
r(0.55yum)/r(350yum) are 1760 and 2080 for MWD and OH, 
respectively. Detailed studies of the possible variations of the 
dust properties are deferred to a later paper. Here the two mod- 
els are used only to test the sensitivity of the results to the actual 
dust properties. The dust emission is calcula t ed with our radia- 
tive transfer program (Juvela & Padoan 2003; Juvela 2005). The 
program uses Monte Carlo simulation to determine the radia- 
tion field intensity at each position within the model cloud. This 



4 The data files describing the dust properties are available at 
http : //www . astro . princeton . edu/ draine/dust/ 
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information is used to determine the distribution of dust tem- 
peratures. The line-of- sight integration of the radiative transfer 
equation then results in synthetic surface brightness maps that 
are calculated for each of the observed wavelengths. 

To avoid edge effects near the boundaries of the area covered 
by the NIR data, we include in the model a larger area than that 
shown in Fig.[T] The modelled area is 10' x 10', centred at coor- 
dinates a 20 oo=19 h O m 54.3 s , £ 2 ooo=-36°55'22". The model will be 
adjusted to reproduce the observed 350 /mi data. This is accom- 
plished by having as free parameters the column densities cor- 
responding to each map pixel. The optimisation will determine 
the column density for each line-of-sight and thus the mass dis- 
tribution in the plane of the sky. In principle, the solution should 
be unique for any combination of the external radiation field and 
dust properties. However, the line-of-sight density distribution 
must be fixed for this modelling. According to Fig.[5]in the plane 
of the sky the FWHM of the column density distribution is ~ 2.5' 
perpendicular to the filamen t. For the assumed distance of 130 pc 
(M arraco & Rydgren 1981) this corresponds to a filament width 
of ~ 0.1 pc. We assume that the filament is approximately cylin- 
drical. Thus, in the model the density distribution along the line- 
of-sight is set to be Gaussian with FWHM equal to 0.09 pc. The 
sensitivity of the results on this assumption is discussed further 
in Sect.ffl The model is modelled by a Cartesian grid of 60 3 cells 
with the cell size corresponding to 10". The surface brightness 
maps are calculated towards one principal axis, the simulated 
map consisting of 60 x 60 pixels. When the model is compared 
with observations, the model results are convolved to the reso- 
lution of the corresponding observation. In particular, at 350yum 
the observations have a resolution of ~ 25 " while the model dis- 
cretisation corresponds to an angular resolution of 10". 

We start by assuming that the filament is illuminated by a n 
isotropic ISRF with the spectrum given by Mat his etaD (1 19831) . 
The solution is obtained iteratively, calculating the model pre- 
diction of the 350//m surface brightness and scaling indepen- 
dently the column density corresponding to each of the 60 x 60 
map pixels. The densities of cells corresponding to a single map 
pixel are scaled with the same number, thus preserving the origi- 
nal shape of the line-of-sight density profile. The first realisation 
is that the CrA cannot be modelled using the standard value of 
the ISRF. The column densities increase without limit but the 
models never reach the surface brightness observed at the centre 
of the filament. There are two solutions to this problem. Either 
the sub-millimetre opacity must be increased by at least factor of 
~2 (or more considering the associated decrease in dust temper- 
atures). The other alternative is to increase the radiation field. 

Figure |7] shows the fit residuals in the case of the OH dust 
model and the scaling of the radiation field by factors ^=2 and 3. 
With^ = 2.0 the surface brightness never reaches the observed 
values, the error at 350 jim remaining above 50%. With^ = 3.0 
the 350 /mi can be fitted but only when the peak Ay values are 
above 100 m over a significant area of the map. Apart from being 
incompatible with the Ay ICER data, the model fails to produce the 
correct shape of the emission spectrum. For example, the 160 yum 
surface brightness is too low typically by more than 20%. 

More satisfactory solutions are found by scaling the ISRF 
values by a factor of x =4. Figure [8] shows the observed sur- 
face brightness maps and the fit residuals for both the MWD and 
OH dust models. Because the model is fitted only based on the 
350yum intensity, and the product of the radiation field and col- 
umn density is capable of producing high enough surface bright- 
ness, the rms error at this wavelength is now below 1%. 

Concentrating on the central area, the fit is better with the OH 
model. Apart from the noise affecting the low column density 
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Fig. 7. The fit residuals for models with the OH dust and the 
radiation field intensity scaled by factors ^=2 and 3. The values 
are relative errors /obs _ ^Mod in per cents of the observed sur- 
face brightness Jobs- The area is this and the subsequent plots of 
model results is 10' x 10' and the plots are centred at coordinates 
^2ooo=19 h O m 54.3 s ,^ 2 ooo=-36°55 , 22". 
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Fig. 8. A model fit to the Herschel surface brightn ess data with 
an is otropic radiation field four times the ISRF of iMathis et al.l 
(1983). The first row shows the original observations. The 
other rows show the errors, (7obs _ ^Mod)A)bs» m un i ts °f P er 
cent. The sec ond row corresponds to the MWD dust model 
(Drain el 120031). and the bott om row to the OH dust model 
(Ossenkopf & Henning 1994). To exclude the regions affected 
by boundary effects, each frame shows the central 8.3' x 8.3' area 
of the whole model that covered an area of 10' x 10 r . To guide 
the eye, the dashed lines show the outline of the area covered by 
the NIR observations. 



part, the relative errors are -10% at both 250 /mi and 500 /mi. 
There is no strong overall bias even at 160yum. This indicates 
that the average shape of the SED is well reproduced. There 
are, however, some systematic effects. The observed level of 
the 500 yum emission is well reproduced in the filament but is 
overestimated is the lower density regions. The error is small, 
-10%, and could be associated with the uncertainty of the inten- 
sity zero point of the observations. A stronger effect is seen at 
the location of the eastern clump, near the boundary of the area 
shown. There the observed 500 jim surface brightness is up to 
-15% lower than the model prediction. The percentage error is 
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not large but it is very clear in the map and significant consider- 
ing that the nearby wavelength of 350 yum is perfectly fitted. In 
the same area the model fails to produce sufficiently high surface 
brightness at short wavelengths, the error increasing to -30% at 
160/mi. One possible interpretation is that the actual radiation 
field is higher in this area. By underestimating the heating, the 
model would require a column density that is too high to explain 
the 350yum map (hence the high 500 yum signal) but does not in- 
clude enough warm dust to reproduce the shorter wavelengths. 
The 160yum (and 250/mi) residuals show a general gradient to- 
wards south-east that can indicate an increase of the radiation 
field intensity. This would be consistent with the higher level 
of NIR scattered light seen in that direction (see Sect. 13.21 ). An 
alternative explanation for the behaviour of the eastern clump 
involves the dust opacity spectral index. A higher ft value might 
be needed to correct the ratio between the 350 yum and 500/mi 
data that are relatively insensitive to variations in the dust tem- 
perature. 

The error maps of the MWD model look qualitatively similar 
to the OH results but the x 2 value (considering maps other than 
the 350//m) is worse by more than a factor of two. Outside the 
filament the 160yum and the 250/mi signals are mainly overes- 
timated and thus the emission appears too warm. However, the 
250-500 yum spectral index of MWD is 2.1, i.e., higher than the 
value y#=1.76 of the OH model.Therefore, compared to the OH 
dust model, the best fit would also be expected to correspond 
to a somewhat lower value of the radiation field. In the east- 
ern clump the 500 yum emission is again overestimated and the 
160yum emission underestimated but, possibly because of the 
different ft values, the errors are smaller. There is still a gradi- 
ent towards south-east consistent with an increase in of the dust 
temperature. However, as already noted, the sign of the 160yum 
errors has changed compared to the OH model, and the MWD 
model would produce insufficient intensity at the southern side 
of the filament. 

We show in Fig. [9] t he Ay map and the filament profile for 
the OH model where the Mat his etaD (Il983l) radiation field was 
scaled by a factor of four. Because of the isotropy of the radi- 
ation field, the shape of the recovered profile is similar to that 
of Ay b_mm . However, the predicted Ay values are not just higher 
than the previous Ay b_mm estimates, which is possible for a num- 
ber of reasons, but also higher than the measurements of Ay ICER . 
This is an indication that either the assumption of the NIR extinc- 
tion curve (used in the conversion of NIR colour excesses to Ay) 
is incorrect or, more likely, the ratio r(0.55 yum)/r(350 yum) used 
in the modelling is incorrect. An increase in the sub-millimetre 
opacity would decrease the model Ay. If this modification were 
restricted to the centre of the filament, the radiation field would 
not have to be changed and the correct shape of the SED would 
still be recovered at low and intermediate column densities. In 
the same figure we also show again the minimum and maximum 
Ay ICER profiles along the selected stripe (dotted lines). The shape 
and magnitude of Ay b_mm in again seen to be in agreement with 
the maximum Ay ICER profile rather than the average profile. 

The previous study of the NIR observations and Spitzer FIR 
data su ggested an anisotr opy with a stronger radiation field in the 
south ( Juvela et aD l2009). Figure [5] lead to a similar conclusion. 
Because the Galactic plane is located north of Corona Australis 
cloud, there must be a more local cause for the asymmetry. The 
star S CrA is located on the south side of the filament, imme- 
diately east of the area covered by our NIR data. The Herbig 
Ae/B e star R Corona Australis is farther in the east dGray et al I 
120031) . The stars appear to affect the dust emission at least in the 
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Fig. 9. The visual extinction in the model with OH dust and an 
isotropic ra diation field with t he intensity four times the stan- 
dard value (M athis et al.lll983l) . The left hand frame shows the 
Ay map from the model smoothed to a resolution of 18". The 
right hand frame shows the Ay profile of the model (dashed line) 
compared to the observed profiles of Ay b_mm (thick solid line) 
and Ay ICER (thin blue line) along the selected stripe. The latter 
two are the same as in Fig. [5j) but are shown at their original res- 
olution. In the plot the Ay 13 " 111111 has been scaled by a factor of 1 .3. 
The dotted lines are the profiles of the maximum and minimum 
Ay ICER along the stripe. 



eastern part of our maps (as seen in Fig. [7]>. The precise spectral 
type of S CrA is not known (the main com ponent S CrA A was 
classified as G5Ve by lCarmona et al.l (120071) ). Based on the ob- 
served K band magnitude of S CrA, — 6.1 m , the NIR inten- 
sity produced by S CrA is comparable to the general ISRF up to 
distances corresponding to several arc minutes. R CrA may have 
similar contribution to the NIR radiation field although, because 
of large intervening extinction, its effect on dust heating is proba- 
bly smaller. The other bright stars in the region, HD 176269 and 
HD 176270, are near the south end of the area covered by our 
NIR observations (see Fig. [2]). The distance estimates of these 
B9V stars are uncertain but they may both be associated with 
the cloud (see lPeterson et al.ll201 ll) . 

We examine the situation by illuminating the model cloud 
with an isotropic ISRF component and adding a second radia- 
tion source in the southern direction, at a position angle of 20 
degrees from south to east. The source is assumed to be at a 
large distance. Both radiation field components are assumed to 
have the same spectral shape as the normal ISRF and their com- 
bined energy input is kept the same as in Fig. four times that 
of the normal ISRF. 

When the southern source stands for 25% of the total radi- 
ation field energy, the north- south gradients in the residuals of 
the 160yum and 500 yum fits disappear (Fig. [TO]). Thus the north- 
south radiation field anisotropy is now correct as far as can be 
concluded from the FIR and sub-mm emission. The result as- 
sumes that the additional radiation is coming directly from the 
SE direction. If the direction for the incoming radiation is not 
perpendicular to the line-of- sight (e.g., more on the front side 
of the cloud), the actual anisotropy of the radiation field may be 
larger. Because the dust is now warmer on the southern side of 
the filament the maximum of the column density should move 
towards north thus increasing the discrepancy with Ay ICER that 
peaks south of Ay b_mm . However, as shown in Fig.[TT] the effect 
on the shape of the model column density profile is not signif- 
icant. The gradient in the east- west direction still remains un- 
explained by the model. More specifically, the clump near the 
east boundary of the NIR map shows, compared to our model, 
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Fig. 10. The relative errors (%) in the predicted surface bright- 
ness in the case of an anisotropic radiation field and the OH dust 
model. 
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Fig. 11. The Ay profile of the filament as determined from 
the radiative transfer model with an anisotropic radiation field 
(dashed lines, scaled with 0.5). The profile is compared with the 

^ NICER (blue soHd lme) and the A sub-mm profile (soHd Wack ^ 

scaled with 1.3). 



an excess at 250 jim and low surface brightness at 500 jim. As al- 
ready noted, this can be a sign of a further temperature gradient 
(natural considering the vicinity of the R CrA region) possibly 
combined with an increase in the dust spectral index. 

If the ISRF level is further increased by 20%, the column 
densities decrease but the x 2 values are higher by more than a 
factor of two, mainly because the 160/mi intensities are overes- 
timated throughout the field. This indicates that, apart from the 
inherent uncertainties involved in the modelling, the total inten- 
sity of the radiation field is well constrained. The models cannot 
produce sufficient surface brightness with lower level of heat- 
ing and with increased ISRF the spectrum becomes incorrect. 
However, if the dust opacity in the cloud centre increases by a 
factor of a few, the lower limit of the ISRF could be relaxed cor- 
respondingly. 

4.2. Modelling of the NIR scattered light 

We calculate the NIR surface brightness maps for the previous 
models where the radiation field corresponded to four times the 
normal ISRF ( Mathi s~et al .11 1 9 8 3l) and the field is either isotropic 
or anisotropic. In the models of the sub-millimetre dust emis- 
sion, the density distributions are similar for the MWD and the 
OH dust, the column densities being only ~ 20% higher for the 
MWD. Therefore, in the following we consider only the den- 
sity distributions derived from the fits of the Herschel data with 
the OH dust model. With this fixed density distribution, the scat- 
tered light can be calculated using both the MWD and the OH 
dusts. The OH model does not specify the scattering function 
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Fig. 12. The scattering and absorption efficiencies for the two 
dust models. The solid curves show the scattering efficiencies 
<2sca for the MWD dust (solid red line) and the OH dust (solid 
blue line) that have been normalised to a value of 1.0 in the V 
band (0.55 /mi). The dotted lines are the corresponding efficien- 
cies for absorption, <2Abs, that have been scaled with the same 
factors as the <2sca values. 



(i.e., the directional distribution of scattered photons) and we 
use the same scattering function as in the case of MWD. 

The differences between the optical and near-infrared prop- 
erties of the two dust models are important as can be seen in 
Fig-US For MWD the scattering is still dominant in the / band 
while for the OH model the efficiencies gsca and gAbs are al- 
most equal. In / band the albedos are 0.68 and 0.58 for MWD 
and OH, respectively. In case of multiple scattering, the effect on 
the surface brightness can be expected to be noticeable. 

The surface brightness images are shown in Fig. [13] For the 
MWD dust, the predicted surface brightness levels are of the 
right order of magnitude. The H band levels are very close to the 
right value, the / band values are -30% higher than the observed 
values and the K band values are lower by a similar amount. This 
could point to a NIR spe ctrum of the ISRF that is re dder than 
in the iMathis et all (Il983h model. lJuvela et al.1 (12008b obtained 
similar results noting that in particular the H and K band values 
were twice as high as expected in the case of the normal ISRF. 
In the present study, the modelling also takes into account the 
effect that the optically thick filament has on the radiation field. 
This shadowing is probably the main reason why the modelling 
prefers higher values for the intensity of the external radiation 
field. 

The NIR calculations fail to reproduce the surface brightness 
dip that is evident in the observations. If the modelling of the 
dust emission underestimates the column density at the centre of 
the filament, this could also partly explain the difference in the 
NIR colours. Because the NIR scattering takes place is mostly in 
the forward direction, an increase in the column density would 
make the spectrum of the scattered light redder. Higher column 
densities would appear to be incompatible with the Ay ICER data 
but, because no stars are seen th rough the centre of the filament 
(see Fig. 7 in I Juvela et al.ll2008l) . this cannot be excluded on the 
basis of Ay ICER alone. The present model is based on the mod- 
elling of dust emission and already has a peak extinction that is 
twice the maximum of the Ay ICER map. If the radiation field in- 
tensity is decreased by a factor of two, the column density could 
increase almost without limit and there would be no constraints 
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on the maximum column density. This is excluded only by the 
fact that a lower intensity of the radiation field would alter the 
SED of the dust emission in a way that is incompatible with the 
observations. 

The OH dust model (Fig. [T3} bottom row) reproduces the 
morphology of the NIR data better than the MWD dust model. 
The dip in the / band is now noticeable at the location of the 
dense clump. The depth of the depression is about half of the 
peak values around the central clump. This is still far from the 
almost complete absence of scattered light in the observations 
but is clearly a step in the right direction. A small depression is 
seen even in the Ks band where the MWD dust produced a clear 
peak. The maximum / band intensity is again slightly higher 
than the observed values. On the other hand, the Ks band signal 
is very low, only one quarter of the observed. 

We show in Fig. [141 similar calculations for the model of 
Fig. [T0| The scattered light is calculated consistently assuming 
that, in addition to an isotropic radiation field component, 25% 
of the total radiation comes from the southern direction. The 
resulting surface brightness asymmetry between the south and 
north side of the filament is already slightly too large. Although 
the ratio of 3:1 between the isotropic and the anisotropic radia- 
tion field components was appropriate for the distribution of dust 
emission, the NIR surface brightness at least precludes larger de- 
gree of anisotropy, assuming that the spectrum of the anisotropic 
source is the same as for the isotropic component. 

The major shortcoming of the models is the absence of a suf- 
ficiently strong surface brightness dip at the location of the main 
clump. This could be caused either by the column densities being 
underestimated (unlikely given the extinction data) or by some 
modification of the dust properties. To check the first alterna- 
tive and to check that the previous differences between the OH 
and MWD models were not just the effect of different opacities 
(instead of the albedo), we recalculated the MWD model after 
scaling its densities by a factor of two. This results only a very 
minor improvement (see Fig. IT3T) although the maximum of the 
Ay map is already well over 100 m . The effect on the morphology 
of the NIR maps is smaller than the difference between the two 
dust models. 

The models predict higher 870 yum intensities than those ob- 
served with LAB OCA (see Fig. [T6l) . However, a large fraction 
of the difference can be explained by the spatial filtering that re- 
moves large scale structure from the LAB OCA map. In the NIR 
observations the / band surface brightness goes almost to zero at 
the centre of the main clump. This is difficult to explain because 
some scattered light is always observed from the outer cloud lay- 
ers. Thus, the explanation may require a very special geometry 
(e.g., a smaller line-of-sight extent of the cloud at the location 
of the clump) and/or further modifications to the dust models as 
suggested by the difference between the MWD and OH model 
results. Although the NIR signal saturates by Ay ~ 20 m , the data 
still set very strong constraints both on the cloud structure and 
the dust properties. 

4.3. Sensitivity to the cloud shape 

In the modelling the column density structure is determined di- 
rectly by the observations. For fixed dust properties and a fixed 
external radiation field the result should be unique. However, 
there are still factors connected with the density distribution that 
can affect the results. The actual three-dimensional shape of the 
cloud is unknown. Although the morphology suggests the pres- 
ence of a single cylindrical filament, the cloud could still be flat- 
tened or elongated along the line-of-sight. Some effects can also 
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Fig. 13. The observed NIR surface brightness (first row) and 
the surface brightness predicted for the models based on the 
Herschel sub-millimetre data (OH dust,^ =4.0). The second and 
the third rows correspond to calculations of the scattered light 
with the MWD and the OH dust models, respectively, in case of 
isotropic illumination. 
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Fig. 14. The NIR surface brightness for the models with 
anisotropic radiation field. The first and the second row corre- 
spond to calculations of the scattered light with the MWD and 
the OH dust models, respectively. 



result from small scale inhomogeneities that allow short wave- 
length radiation to penetrate deeper into the cloud. We tested the 
sensitivity of the results to these factors using the model pre- 
sented in Fig.fTOl 

The dumpiness was implemented by multiplying the den- 
sity of each cell in the model by a random number uniformly 
distributed between zero and one. This results in a fair amount 
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Fig. 15. The observed NIR surface brightness (first row) and 
the surface brightness predicted by models with the densities 
twice the values of those fitting the sub-millimetre data. The sec- 
ond and the third row correspond to the cases of isotropic and 
anisotropic radiation field, respectively. 
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Fig. 16. Comparison of the observed 870 /mi surface brightness 
and the prediction from the model with x - 4.0 with the OH 
dust. 



of inhomogeneity considering that the cloud radius is only -30 
cells. All densities were then rescaled so that the model again re- 
produced the observed 350 yum surface brightness. Compared to 
the original model with a smooth density distribution, this re- 
sulted in less than 10% changes in the column density along 
individual lines-of-sight. More importantly, this level of inho- 
mogeneity did not introduce significant systematic change in the 
column density nor in the intensity of NIR scattered light pro- 
duced. Of course, the effects would be stronger if the cloud in- 
cluded large continuous cavities extending to the centre of the 
cloud. 

The overall cloud shape was varied by changing the FWHM 
of the line-of-sight density distribution by ±30%. The column 



densities were again adjusted according to the 350 /mi data. 
When the cloud was 30% more extended along the line-of-sight, 
the resulting column density towards the model centre was lower 
by -20%. When the cloud was flattened by 30%, the percentual 
increase of the model column density was closer to 40%. These 
numbers do not yet take into account the fact that a for the best 
fit to the observed SED also the intensity of the external radia- 
tion field might have to be readjusted. This could reduce the net 
change of the column density. However, with the radiation field 
of Fig.[lOlthe 160yum surface brightness was already correct to 
within -10%. 

We conclude that the uncertainty of the three-dimensional 
shape of the cloud translates to an uncertainty of some tens of 
per cent in the column density. These uncertainties affect the ac- 
curacy to which the dust opacity can be determined. However, 
the modelling is better constrained when independent extinction 
measurements are available. 

5. Conclusions 

We have examined the structure of the northern filament of the 
Corona Australis cloud using the combination of Herschel sub- 
millimetre data and near-infrared observations. The study has 
lead to the following conclusions: 

- At Ay < 15 m the three column density estimators, the ther- 
mal dust emission, the NIR reddening of the background 
stars, and the surface brightness caused by the NIR scattered 
light, are all linearly related to each other. The differences 
in the actual numbers can be explained by the uncertainty of 
the dust opacity and of the radiation field intensity. 

- At low column densities (Ay ~ 10 m or below) the scattered 
light provides the best column density map in terms of reso- 
lution. Its predictions are very tightly correlated with the es- 
timates obtained from the dust emission. However, the slope 
of the relation depends on the intensity of the radiation field. 

- Based on the dust emission, the north-south column density 
profile was observed to be skewed with a sharp drop in the 
column density on the northern side. The profile obtained 
from the reddening of the background stars is more symmet- 
ric but is strongly affected by the lack of stars visible through 
the central filament. 

- The modelling of the dust emission suggests that the radi- 
ation field intensity at the location of the cloud filament is 
about four times the value of the normal ISRF. The value 
can be lowered only by assuming that the dust sub-millimetre 
opacity has increased at the centre of the filament. However, 
the models of the NIR scattered light are also consistent with 
an elevated radiation field intensity. 

- According to the models, the radiation field is anisotropic 
with an approximate ratio of 3:1 between the isotropic com- 
ponent and additional radiation coming from the southern 
direction. 

- The models are unable to reproduce the deep dip in the NIR 
surface brightness at the centre of the filament. One possible 
explanation is a change in the dust properties (as indicated 
by the differences between the dust models examined) that 
has lowered the NIR albedo. 

Acknowledgements. MJ acknowledges the support of the Academy of Finland 
Grants No. 127015 and 250741. 
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